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Abstract 

The complete elastic tensor of mullite has been determined by brillouin spectroscopy at 
room temperature and elevated temperatures up to 1200 °C. Equivalent, isotropic moduli (bulk, 
shear, and Young’s) have been calculated. The room temperature values obtained using Voigt- 
Reuss-Hill averaging are: Kvrh = 173.5J- 6.9 GPa, Gvrh ~ 88.0 ± 3.5 GPa, E V rh = 225.9 ± 9.0 
GPa. All moduli show relatively gradual decreases with temperature. The temperature 
derivatives obtained for the equivalent, isotropoic moduli are: dKyRH^T = -17.5 + 2.5 MPa/ C, 
dG VRH /dT = -8.8 ± 1 .4 MPa/°C, dE VRH /dT = -22.6 ± 2.8 MPa/°C. Substantial differences 
between bulk properties calculated from the single crystal measurements in this study and the 
properties reported in the literature for polycrystalline sintered mullite are identified, indicating 
the importance of factors such as microstructure, intergranular phases, and composition to the 
elasticity of mullite ceramics. 


” - Fellow American Ceramic Society 


1 



> 


I. Introduction 

Mullite is found in many common ceramics. More recently it has gained favor for use in 
advanced materials such as structural ceramic matrix composites as both a matrix and reinforcing 
phase. It has such desirable qualities as high temperature stability and oxidation resistance that 
give it potential for use in many applications 1 such as gas turbine engines . 2 

Because of the long-term use of mullite and its technological importance , 3 many 
measurements have been made of its mechanical properties, including measurements of its elastic 
properties at ambient and elevated temperatures . 4-10 These have all concentrated, however, on 
bulk polycrystalline properties, and there are great discrepancies in the measured values. 

Recently, mullite has become available in single-crystal or highly textured polycrystalline 
fiber form by use of a laser-heated float zone crystal growth technique . 1 1 Such samples allow 
measurement of the single-crystal elastic properties of mullite. These measurements provide 
useful property data for the design of components using single crystals (e.g. single crystal 
reinforcing fibers). They also yield information as to what degree the elastic properties of mullite 
ceramics are a result of the intrinsic properties of the mullite crystal as compared to 
characteristics that may depend on processing such as microstructure, and intergranular phases 
which are common in sintered mullite. Therefore, a detailed understanding of the single ciystal 
elastic properties of mullite is extremely useful for a fundamental understanding of the mechanical 
properties of complex, alumino-silicate ceramics. Such single crystal measurements are made in 
this study using Brillioun scattering. 

Recently, we have reported a full single crystal elastic tensor for mullite at room 
temperature and selected constants (Cj b C33, C 44 , C 55 ) to 1400°C using brillouin spectroscopy . 12 
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Orientation of the samples in the previous work was based on sample morphology and 
symmetry of the velocity data itself, and the previous study utilized a large rigidly mounted 
furnace that allowed access to only a few phonon directions with limited knowledge of the 
sample orientation at high temperature. This work gives a more accurate room temperature 
tensor than the previous work, due to a refinement of the sample crystal orientations by X-ray 
diffraction, as well as the temperature dependence of the full tensor to 1200 °C also with 
improved accuracy made possible by the use of a compact furnace design that allows collection 
of multiple phonon directions without remounting the sample. 13 Neutron diffraction to 1650 °C 
and X-ray diffraction to 900 °C of polycrystalline or unconstrained powdered mullite with a 
composition of 3 Al 2 0 3 -2Si0 2 was also reported in the previous work. This data was used to 
calculate thermal expansion for the previous study. For the current study, thermal expansion 
was taken from an X-ray study of mullite with composition 2 Al 203 -SiC >2 to 900 °C because of its 
closer composition to the current sample. 14 

II. Theory 

Brillouin scattering is the inelastic scattering of light from acoustic waves, allowing the 
determination of the wave speed. It provides an ideal probe of the elastic properties of the single 
crystal mullite samples available (especially at high temperatures), requiring only small sample 
volumes (e.g. 100 x 100 x 50 pm) and no physical contact with the sample. It is capable of 
measuring the phonon velocities along many propagation directions in a single sample. Brillouin 
scattered light is shifted in frequency with respect to the incident light by a factor that is 
proportional to the velocity of the acoustic waves (Eq. I). 15 
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( 1 ) 



Here V is the velocity of the acoustic wave, Aco is the frequency shift of the scattered light, co is 
the frequency of the incident light, c is the speed of light, n is the index of refraction of the 
sample, and 0 is the scattering angle. 

For this study, platelet geometry is used such that the sample has flat, parallel faces, the 
phonon propagates in the plane of the sample, and the face normals, incident, and scattered light 
rays all lie within a plane. 16 Figure 1 gives a schematic of the platelet scattering geometry. With 
the symmetric platelet geometry, the scattering angle is easily defined, and no knowledge of the 
index of refraction of the sample is necessary. This is particularly beneficial for high temperature 
studies where the variation of index of refraction with temperature is generally unknown and 
difficult to measure. 

The elastic waves probed may be treated by a continuum mechanics analysis that relates 
wave velocities along a given crystallographic direction to the elastic moduli and density of the 
material via Christoffel’s equation (Eq. 2). 17 
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C is the elasticity tensor for the material, q is a unit vector along the path of phonon propagation, 
p is the density, and 6 is the Kronecker delta function. Since mullite is orthorhombic, 9 
independent elastic constants (Cn, C 22 , C 33 , C 44 , C 55 , C^, C[ 2 , C 13 , C 23 ) are required to describe 
its elastic tensor (represented here in Voigt notation). 


III. Experimental Procedure 

(1) Brillouin Scattering System 

For this study, an Ar + laser with a wavelength of 514.5 nm is used to illuminate the 
sample. Scattered light is analyzed by a 6 -pass tandem Fabry-Perot interferometer. Figure 2 is a 
schematic of the system used to collect Brillouin spectra. The operation of the interferometer 
has been detailed elsewhere. 18 - '9 

Room temperature work is performed with a 90° scattering angle. An Eularian cradle is 
used to rotate the sample around its face normal to access different phonon directions. The high 
temperature results of this study are obtained using a compact furnace mounted on the Eularian 
cradle, allowing multiple phonon directions to be collected without remounting the sample. The 
construction and operation of this furnace is described elsewhere . 13 An 80° scattering angle is 
used for the high temperature work. 

(2) Sample Preparation 

Samples are taken from a fiber grown using the laser-heated, float zone method 11 . For 
preparation of the source rod, high purity (99.99% pure) polycrystalline alumina powder was 
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obtained from CERA C/pure (325 mesh) and 99.99 % pure SiC >2 from Alpha Products. The 
starting composition was chosen close to the 3:2 molar ratio of AI 2 O 3 : SiC> 2 and was 70.4 wt % 

AI 2 O 3 . The fibers characterized were from a second pass through the laser. Figure 3 shows side 

and cross sectional views of a representative mullite fiber. A definite growth face is apparent, 

forming a sawtooth pattern. X-ray analysis shows this to be the (010) face. Cleavage is 

occasionally observed along this plane, as is observed in natural mullite and sillimanite samples . 20 

The fiber has a roughly elliptical cross section minus the growth facet. The roughness on the 

surface in Figure 3 (b) is due to the diamond saw cutting process. The fiber is for the most part 

optically clear, with occasional white streaks across its width, apparently consisting of bubbles 

or occlusions. These white streaks often coincide with the points of the sawtooth pattern along 

the fiber edge. 

Separate platelets for velocity measurements with approximate face normals of ~[0 1 0], 
~[0 01], and ~[0.9744 0.2250 0] were ground using a series of diamond and SiC abrasives ranging 
from 30 to 1 pm. (Normals and planes are given in an equidistant Cartesian coordinate system, 
not the orthorhombic lattice of mullite.) Machine grinding was used for some of the coarser steps 
and the finer polishing was completed by hand. Glass was placed around the sample as a 
sacrificial material to minimize rounding of the sample faces, which was less than 1.5° for the 
entire included angle of the face (usually < 1 °). Nonparallelness between the faces was less than 
0.5°, and the faces produced had a mirror finish with very few scratches visible at 40x 
magnification. Final dimensions of the samples taken lengthwise from the fiber, ~(100) and 
-(010), were approximately 3 mm x 0.4 mm x 0.25 - 0.1 mm. Cross sectional samples, ~(001), 

* - Ceralox Corp., Tucson, AZ 08576 
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had approximate diameters of 0.4 mm and thicknesses of 0.15 -0.1 mm. The orientations of the 
sample platelets were determined using X-ray precession photography. 21 


IV. Results 

(!) Sample Characterization 

The fiber shows uniform extinction under cross polarized light when the polarization 
direction is nearly aligned or perpendicular to the fiber axis, regardless of rotation about this axis. 
When ground with parallel faces, some variation in interference color is noted occurring in bands 
across some samples. This may be a result of a slight change in index of refraction due to 
compositional variations (or possibly slight variations in thickness). Extinction is somewhat less 
uniform in cross sections of the fiber. 

The composition of the mullite was determined with an electron microprobe (wavelength 
dispersive) using natural kyanite as a standard. The average composition is listed in Table I along 
with unit cell parameters determined by X-ray diffraction for the mullite sample. The error 
quoted for composition is the standard deviation of values measured at 10 pm intervals along the 
length of a fiber section and 5 pm intervals across its width. Density is determined using Eq. 3 
for the chemical formula of the unit cell. 22 


Al2(Al2+2 x Si2-2x)0 1 0 -x 


( 3 ) 
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x is the average number of oxygen vacancies per unit cell. It is related to the Al 2 0 3 /Si02 molar 


ratio, n(Al 2 0 3 /Si0 2 ) by Eq. 4. 


2n(Al 2 0 3 /Si0 2 )-2 
2n(Al 2 0 3 /Si0 2 ) + 1 


Initial X-ray diffraction analysis was performed using a four-circle diffractometer and 
shows that the sample is actually a very highly textured polycrystal, as seen in Figure 4. This 
analysis verifies that the crystallites have orthorhombic symmetry with slightly varying 
orientations defined by small rotations about the c axis, which is coincident with the fiber axis. 
All the crystallites found lie within 5° of each other and the majority lie within 2°. There is no 
resolvable splitting in the c axis peaks. Lattice parameters for the sample are calculated by fitting 
the major peak for each scan. A total of six peaks are used in the calculation: (040), (400), 

(002), (402), (042), (222) Jhe errors quoted encompass all cell parameters calculated with five 
of the six peaks (i.e. the maximum and minimum values from six calculations with five peaks). 

The calculated unit cell angles deviate slightly from 90° (a = 90.4174°, |3 = 90.2971°, y = 
89.8230°), probably due to the small peaks in close proximity resulting from the highly oriented 
polycrystalline nature of the sample and large sample size. For the purpose of calculating the cell 
volume, the interaxial angles were fixed to a = p = x = 90°. The sample characterized in Table I 
is the same as that in the previous study, but the inclusion of some additional X-ray data resulted 
in a slight (0.6%) decrease in calculated volume. 

(2) Room Temperature Elasticity 
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The full set of 9 room temperature, elastic moduli (Table II) have been determined using a 
linearized least- squares inversion procedure based on 158 separate velocities measured along 37 
different crystallographic directions in three different planes(Figures 5, 6, and 7). This algorithm 
calculates velocities based on an intial set of preliminary Cy s (Eq. 2) along the same phonon 
propagation directions as those measured and then varies the elastic moduli so as to minimize the 
misfit between the calculated and measured velocities.^ V elocities calculated from the moduli 
obtained appear as solid curves in Figures 5, 6, and 7. We note that the final best-fit solution is 
highly insensitive to the starting guess for the Cy’s. Equivalent, isotropic moduli have been 
calculated using the Voigt, Reuss, and Voigt-Reuss-Hill averaging schemes (denoted by V, R, and 
VRH subscripts) and are listed in Table II. 

The errors in Table II refer purely to the statistical error in the fit. Deviations from the 
nominal sample normals (~[0 1 0], ~[0 0 1], -[0.9744 0.2250 0]) are less than -4° for all samples. 
The higher symmetry sample normals are retained to simplify calculation of moduli and error 
analysis. Rotation about the sample normal is determined by extinction under cross polarized 
light and is considered accurate to +3°. The error resulting from misorientation of the sample 
(including rotation about the sample plane normal) is estimated to be less than + 2% of the 
calculated moduli. This estimate is based on including deviations in each axis of the sample plane 
normal and each rotation singly and assuming additivity. (Rotations are calculated with sample 
plane normals of [1 0 0], [0 1 0], and [0 0 1] for simplicity.) The uncertainty in density, on 
which the moduli are directly proportional, is ± 0.7%. 

(3) High Temperature Elasticity 
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For high temperature moduli determinations, velocities were collected from each sample 
plane at 200 °C increments to 1200 °C. Figures 8, 9, and 10 show the calculated elastic 
longitudinal (e.g. C| i), shear (e.g. C 44 ), and cross (e.g. C 12 ) moduli along with best fit, least 
squares regression lines. High temperature constants are calculated at each temperature from 16 
velocities collected along 9 directions (3 in each sample plane; one along each projection of the 
principle axes and one 45° between them). Linear temperature dependencies are suggested for all 
single crystal moduli, with the possible exception of C 2 3 which shows a large scatter with 
temperature. Linear decreases in E V r», Gvrh> and K V rh are also indicated (Figures 11, 12, 13). 

The high temperature constants are calculated directly from Eq. 2 using the same sample 
face normals as the room temperature calculations. C 12 and C 13 are determined from a weighted 
average of a single shear and longitudinal acoustic modes based on the sensitivity of the moduli to 
the specific velocity. 21 C 23 is calculated only from the longitudinal mode due to difficulty in 
assigning a specific shear mode. This may account for some of the scatter in the data. A volume 
coefficient of thermal expansion of 15.8 x 10' 6 /°C for 2A1 2 0 3 -Si0 2 mullite is used for the 
calculation of density at high temperature as obtained from the sum of linear expansion 
coefficients of the lattice parameters given by Schneider and Eberhard. 14 This value was 
determined from measurements only to 900 °C. Because of the lower temperature limit for this 
study (1200 °C) as compared to the previous one (1400 °C), it was decided that the extrapolation 
in temperature is more acceptable than the greater difference in composition. The temperature 
derivatives of the single crystal moduli and Voigt- Reuss-Hill average moduli are given in Table III. 
Errors quoted refer to the 95% confidence intervals for the slope of the linear fit. 
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V. Discussion 


(1) Room Temperature Elasticity 

The single-crystal moduli, for the most part, fall within the error bars of those previously 
reported. 12 The one significant discrepancy is C 23 which was previously reported as 135 GPa 
based an orientation resulting from the symmetry of the velocity data itself. The current 
reorientation represents a rotation of 13° (around the c axis) for one plane, ~(100), which strongly 
affects C 23 . The current value (142.3 GPa) is ~5% higher and is considered the more accurate due 
to the precision of the X-ray orientation. The other constants are clearly affected much less by 
this reorientation. This occurs since the other constants depend very weakly on this rotation 
(Eq. 2). Cross constants such as C 23 are generally the most strongly affected by 
misorientation. 24 

Sillimanite, Al 2 Si0 5 , which has a very similar crystal structure to mullite 25 , shows similar 
elastic behavior. The relative magnitudes of the moduli by value is the same for both crystals, 
except for the exchange of Ci 3 and C'gg. Andalusite shows the same order in longitudinal 
constants when its a and b axes are exchanged, which is the transformation required to bring its 
structure into closest correspondence with mullite 25 . The bulk properties calculated for mullite 
and sillimanite (Table II) are similar, while andalusite has a smaller bulk modulus but larger shear 
and Young’s moduli. 

Bulk values calculated here for mullite agree well with the recent measurements of a 
polycrystalline mullite sample by Ledbetter et al. 19 , as well as the polycrystalline Young s 
modulus measurements of Pentry 2 . There are, however, considerable differences among some of 
the values of elastic moduli reported for polycrystalline mullite (Table IV), and certain authors 8>9 
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have reported values that exceed even the Voigt bounds of the respective moduli calculated here. 
These should set an upper limit on the possible values for a polycrystalline aggregate of mullite. 
The discrepancy may be due to intergranular phases. Compositional differences may also play a 
role. The polycrystalline studies were performed on mullite samples with a lower alumina 
content (3 A1 2 0 3 • 2 Si0 2 ). 

(2) High Temperature Elasticity 

Comparisons with temperature derivatives of polycrystalline elastic properties given in 
the literature show significant differences (Table V). The derivative obtained for Young’s 
modulus in this study is more negative than -8 MPa/°C as seen in the data of Wachtman for 
temperatures below 600°C 5 . (Note that the scale in Figure 1 1 for the literature data is twice as 
fine than for the current data, while in Figure 12 it is twice as coarse.) Ledbetter also reports a 
more gradual decrease in Young’s modulus (-11 MPa/°C) 10 . These differences could result from a 
glassy phase in the polycrystalline body. Silica glass is known to stiffen with temperature 26 . A 
similar effect may occur with an intergranular glass phase, thereby mitigating the reduction in 
stiffness of the mullite grains. The temperature derivative of the shear modulus from the present 
study, however, is much smaller than -30 MPa/°C which is suggested by a linear fit to Wygant’s 
data. 4 This may be due to a low melting (or more rapidly softening) intergranular phase. 

Ledbetter reports an extremely gradual -3.3 MPa/°C temperature derivative of shear modulus as 
well as a bulk modulus temperature derivative of -20 MPa/°C 1 ®. All values from Ledbetter are 
for room temperature and are corrected to zero porosity. The temperature behavior of mullite 
ceramics is clearly affected, and possibly dominated, by properties other than those intrinsic to 
the crystalline mullite itself. 
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VI. Conclusions 


The full elastic tensor of mullite has been determined at room and elevated temperatures 
to 1200 °C. Room temperature values show close correspondence with structural analogs, 
sillimanite and andalusite. Bulk elastic properties calculated from the single crystal moduli agree 
well with some measurements on polycrystalline mullite 7 ’ 10 , but there are substantial differences 
with much published data. In particular some authors report values 8 - 9 that exceed the maximum 
possible limit (Voigt bound) for average moduli calculated here, suggesting that intergranular 
phases and composition may play a very large role in determining the elastic properties of 
nominally pure mullite ceramics. 

High temperature measurements were made in 200 °C intervals to 1 200°C. A relatively 
gradual linear decrease was noted in all moduli with temperature. The temperature derivatives of 
bulk properties calculated from the single crystal moduli do not correspond well with much of 
the data in literature, again suggesting that intergranular phases are important in the elasticity of 
mullite ceramics. Phases having behavior at high temperature much different than mullite, such as 
stiffening or rapid softening with temperature, are suggested by this data. 
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Figure 1 : Schematic of platelet scattering geometry 
Figure 2: Schematic of Brillouin scattering system. 

Figure 3: Optical micrographs of mullite sample (a) lengthwise and (b) in cross section 

Figure 4: X-ray scan rotating sample around fiber axis with fixed 20 

Figure 5: Room temperature velocities versus rotation in ~(100) plane 

Figure 6: Room temperature velocities versus rotation in ~(010) plane 

Figure 7: Room temperature velocities versus rotation in ~(001) plane 

Figure 8: Single-crystal longitudinal moduli versus temperature 

Figure 9: Single-crystal shear moduli versus temperature 

Figure 10: Single-crystal cross moduli versus temperature 

Figure 1 1 : Young’s modulus of mullite versus temperature 

Figure 12: Shear modulus of mullite versus temperature 

Figure 13: Bulk modulus of mullite versus temperature 
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Table I: Properties of mullite sample at room temperature 


a (nm) 

0.758 + 0.001 

b (nm) 

0.767 ±0.001 

c (nm) 

0.289 ±0.001 

Cell Volume (nm 3 ) 

0.168 ±0.001 

Molar AI 2 0j/Si02 ratio 

2.5 ±0.1 

Theoretical Density (g/cm 3 ) 

3.11 ±0.02 
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Table II: Room temperature moduli of mullite, andalusite, and sillimanite 


Constant 

Mullite 

C„ (GPa) 

281.9 ±0.7 

C 22 (GPa) 

244.2 ± 0.9 

C 33 (GPa) 

363.6 ± 1.0 

C 44 (GPa) 

1 11.7 ± 0.6 

C 55 (GPa) 

78.2 ±0.5 

C 66 (GPa) 

79.2 ± 0.5 

C 12 (GPa) 

105.1 ± 1.0 

C13 (GPa) 

100.3 ± 1.3 

C 23 (GPa) 

142.3 ± 1.6 

K v (GPa) 

176.1 

K R (GPa) 

171.0 

Kvrh (GPa) 

173.5 

G v (GPa) 

90.0 

Gr (GPa) 

86.1 

Gvrh (GPa) 

88.0 

E v (GPa) 

230.6 

E r (GPa) 

221.2 

Evrh (GPa) 

225.9 


Sillimanite 27 Andalusite 27 


287.3 

233.4 

231.9 

289.0 

388.4 

380.1 

122.4 

99.5 

80.7 

87.8 

89.3 

112.3 

94.7 

81.4 

83.4 

116.2 

158.6 

97.7 

175.1 

165.9 

166.4 

158.0 

170.8 

162.0 

95.1 

100.4 

87.8 

97.7 

91.5 

99.1 

241.6 

250.6 

224.0 

243.0 

232.8 

246.8 
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Table III: Fitted temperature derivatives for mullite to 1200°C 


Modulus 

Temperature derivative (MPa/°C) 

Cn 

-29.2 ±3.2 

C 2 2 

-29.7 ± 1.6 

^33 

-28.6 ±4.2 

C44 

-7.9 ± 1.1 

C55 

-6.5 ±0.9 

C66 

-10.1 ±0.9 

Cl2 

-16.1 ±0.5 

C[3 

-6.8 ± 1.4 

^23 

-8.1 ± 10.5 

Kvrh 

-17.5 ±2.5 

Gvrh 

-8.8 ± 1.4 

Evrh 

-22.6 ±2.8 
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Table IV: Bulk properties of polycrystalline mullite reported in literature 


E (GPa) 

V 

G (GPa) 

K (GPa) 

Density 

Ref. 

225.9 


88.0 

173.5 

3110 kg/m 3 

Pres. Study (VRH 






Aver.) 



59 



4 

143.0 




2780 kg/m 3 

5 

124 





6 

223.4 ± 5.5 



Corrected to 1 00% 

7 





(Spriggs) 


224.1 ±5.5 



Corrected to 100% 

a 





(Hasselman) 


246 




>98% 

8 





(3150 kg/m 3 ) 


194 

0.276 

76.0 

144.3 

95% 

28 

252 

0.266 

99.5 

179.5 


9 

227.5 

0.280 

88.9 

172.4 

99.6% 

10 

229.1 

0.280 

89.5 

173.9 

Corrected to 100% 

a 
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Table V: Comparison of temperature derivatives with literature values 


dK/dT (MPa/°C) 

dG/dT(MPa/°C) 

dE/dT (MPa/°Q 

Ref. 

Notes 

-17.5 ±2.5 

-8.8 ± 1.4 

-22.6 ±2.8 

Pres. Study 

VRH average 





20- 1200 °C 



-8 

5 

20 - 600 °C 

-20 

-3.3 

-11 

10 



-30 


4 
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Figure 2: Schematic of the Brillouin spectroscopy system (BS = beam splitter, M = mirror, L = 
lens, FL = focusing lens, CL = collecting lens, PH = pinhole, PR = prism, FPC = Fabry-Perot 
controller, OSC = oscilloscope, MCA = multichannel analyzer, PC = Personal Computer, FPP = 

Fabry Perot Plates) 
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1800 



Figure 4: X-ray scan rotating 
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Figure 5: Room temperature velocities versus rotation in ~(100) plane 
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Figure 7: Room temperature velocities versus rotation in ~(001) plane 



400 


o 

o 

o 

o 

o 

o 

O 

O 

O 

O 

CO 

0> 


<N 

o 

CO 

vo 


<N 

O 

cn 

cn 

cn 

cn 

cn 

<N 

(N 

CN 

(N 

<N 


(bjo) sn[npoj\[ 





(EJO) siqnpoi\[ 


Figure 9: Single-crystal shear moduli versus temperature 
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Figure 11: Young’s modulus of mullite versus temperature 
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Figure 12: Shear modulus of mullite versus temperature 
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Figure 13: Bulk modulus of mullite versus temperature 




